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Abstract 

To obtain the value of electromagnetic coupling constant at q 2 = m z , a, which 
plays a key role in electroweak physics one has to integrate the cross-section of e + e~- 
annihilation into hadrons divided by (s — rn 2 z ) over s from threshold to infinity. By 
combining, for each flavor channel, the contribution of lowest resonance with the 
perturbative QCD continuum, we obtain 1/a = 128.90 ± 0.06 a result which is close 
to known result obtained with purely experimental inputs, i.e. 1/a = 128.87 ± 0.12. 
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The detailed analysis of electroweak observables starts from three input parameters: G^, 
the Fermi coupling constant (extracted from muon decay), mz, the Z-boson mass (measured 
at LEP) and a, the electromagnetic coupling constant at q 2 = m z , obtained from dispersion 
relations. In fact, a Born aproximation to the minimal standard model which starts with 
a (rather than a = a(0) = 1/137.0359895(61) ) reproduces the precise experimental values 
of the Z-decay parameters (obtained at LEP) and of the W mass (obtained at hadron 
colliders) with unexpectedly high accuracy [|TJ, 0. For example for the ratio of vector 
and axial coupling constants of the Z-boson to charged leptons one obtains in this at Born 
approximation 0: 



[gv/g A \a = 0.0753(12) , (1) 
while the latest experimental numbers are ||: 

\9v/9a\lep = 0.0711(20) , (2) 

\9v/9a\lep+sld = 0.0737(18) . (3) 
If instead of a one uses a(0), then one gets: 

[gvlg A ]a = 0.152 , 

which is about 40a's away from experiment as was stressed in 0. The value of a is of fun- 
damental importance, and its error determines the uncertainty in the theoretical prediction 

©• 

a is defined through the following formulas: 

5a = 2,(0) - , (5) 

' m z 

where in @ charged leptons and five quark flavor contributions in photon polarization 
operator should be taken into account. Contributions of (if) and (VyVF) loops may be 



omitted in (|5|); they are numerically small and usually are attributed to proper electroweak 
radiative corrections H. The following integral representation for 5a is valid: 

5a = / ^+e r an» ds } (6) 
4ir z a J m z — s 

where integral goes from threshold to infinity and its principal value at s — m\ should be 
taken. The lepton contribution of e, /i, and r to (||) are readily calculated and one gets: 

(5a)! = ^[S In ^| - | = ^[22.5 + 11.8 + 6.2] = 0.0314 . (7) 
Sir mj 6 Sir 

For the hadronic contribution in j5| the following number was obtained (see also ||): 

(5a) h = 0.0282(9) . (8) 

To obtain this number the experimental cross-section for e + e~-annihilation into hadrons 
below s = (40GeV) 2 and parton model result above s was used in and 0. 

The difficulty in the theoretical determination of (5a) h comes from its logarithmic de- 
pendence on the infrared cutoff. As it was mentioned in 0, the result of the dispersion 
calculation of (5a)h can be reproduced by using perturbative QCD with the following effec- 
tive " quark masses" : 

m u = 53MeV , m d = 71MeV , m s = 174MeV , 

m c = 1.5GeV , m b = 4.5GeV . (9) 

Unfortunately one can not attribute any physical meaning to these values of m u and md- 

Our aim here is to present simplest sensible model for cr e + e -_> hadrons which can simulate 
the result given in (|8]). To do this we use one physical resonance (p,u,(p, J/\l/ and T) at 
the beginning of spectrum and then starting from = rrii + ^ the QCD improved parton 
model continuum in each quark channel. 

For resonance contribution we use Breit-Wigner formula: 



Substituting in (|6|), neglecting terms of the order of (m/rriz) 2 and integrating from — oo to 



m + ^ we obtain: 



iuOL ) resonance * • (H) 

a m 4 



Thus vector meson contributions into 5 a are: 

H * 1 (12) 

5a 0.00274(13) 0.00024 0.00042 0.00053 0.000045 , 

where we take into account experimental uncertainty for p-meson contribution as the only 
noticeable. 

For continuum contribution we use the following formulas: 

^=l = 27T^(l + ^) , (13) 
S 7T 

2n a 2 ,„ aJs). 

= ?7* + ir> • (14) 
= w (1 + o^Oj (15) 



167r« 2 / 4m 2 2mJ w , a s (s) x 
a ce = Wl ^ 1 + — * 1 + , 16 

9 S V S S 7T 
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where we use for a s (s) the following formula: 

12vr 

« s (s) 



(33 -2n/)lns/A( n /) 2 v 7 

with a s (mz) = 0.129(5) as an input (this one loop value corresponds to 0.125(5) at three 
loops, which is extracted from latest LEP data 0). We take rif = 5 for s > my, rif — 4 for 
m\ > s > m 2 ^, rif = 3 for m 2 ,^ > s > m 2 and rif = 2 for m 2 > s > (m p + T p /2) 2 . This 
corresponds to A (5) = 160MeV, A< 4 > = 220MeV, A (3) = 270MeV and A (2) = 300MeV. 



Substituting (|T3|) - ( |18D into (Jg) with m c = mj = we get: 



1=1 1 = ss cc bb 

(19) 

5a 0.01174 0.00133 0.00249 0.00741 0.00123 
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Summing up contributions of ([12]) and fll9D we get: 

(5a) h = 0.0282 , a = (128.87)" 1 (20) 

Comparing with obtained by integrating experimental data results (|]) (5a)h = 0.0282 and 
a = [128.87(12)] _1 we see that agreement is astonishing. The contribution of a s correction 
in ([19]) is rather small, 0.00087 + 0.00010 + 0.00018 + 0.00042 + 0.00006 = 0.00163, so even 
if light gluino octet slow down a s running in order to accomodate a s values measured at 
quarkonium decays || (5a)h will decrease by 0.0002 only. 

We have to make a few comments: 

(1) taking contributions ~ a 2 s in continuum cross-section into account and using next-to- 
leading order formula for a s (s) we increase (5a)h by 0.00045; negative contribution of ~ a^ 
term appeares to be approximately two times larger. In beauty and charm channels third 
loop gives much smaller contribution than second (numerically both are negligible) so we 
can trust our continuum calculation. In strange channel third loop contribution equals that 
of second, while in / = 1 and 1 = channels it is two times larger. So below, say, 1.5GeV 
perturbative continuum can not be approved. Allowing physical continuum variation at the 
level of ±15% around tree plus one loop perturbative continuum value in the domain 1 - 
2GeV we get ±0.0004 variation in (5a) h', 

(2) experimental uncertainty in Tu of vector resonances lead to (5a)h variation of the 
order of 0.0002, while that in a s (mz) - to 0.0001 variation. Both are small compared with 
the uncertainty 0.0009 in @; 

(3) subtracting from the p contribution the integral over Breit-Wigner formula from — oo 
to two pion threshold we diminish it by: 

3r 

5a sub = —2 arctan TJ(2(m p - 2m*)) = 0.00017; (21) 

27iam p 

(4) taking into account heavy quark masses m c = 1.6 GeV, m;, = 4.7 GeV, we decrease 
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(Sa)h correspondingly by: 



(5a h ) m = 0.00031 + 0.00008 = 0.00039; (22) 

(5) finally, at energies E = + y our model curve for cr e + e - ^hadrons is discontinuous. 
To understand (Sa)h sensitivity for the details of the model we change it in the following 
way: we continue p, u, tp and J/\I/ resonance curves up to their intersection with quarks 
continuum. In this way [5a)h increases : 

S(5a) h = 0.00039 (23) 

Subtracting from (|23| ) sum of fl22|) and ([H]) and taking uncertainty from point (1) above 
for total shift we get: 

(5a) h = 0.0280(4), a = (128.90(e))" 1 (24) 

So it is evident that the value of (5az)h is rather insensitive to the details of the model of 
<7 e + e -^hadrons- More refined model which takes all known resonances in each flavour channel 
into account gives (Sa)h = 0.0275(2) |J. 

For real progress in diminishing error in (§) systematic error in cross section of e + e~- 
annihilation into hadrons in background region below 3 GeV should be improved || , [|Kj . 
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